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Wc present measurements of branching fractions for five B-mcson decays to two-body charmless 
final states. The data, collected with the BABAR detector at the Stanford Linear Accelerator Center, 
represent 459 million BB pairs. The results for branching fractions are, in units of 10~® (upper 
limits at 90% C.L.): B{B+ ?7p+) = 9.9 ± 1.2 ± 0.8, B{B° -> rj'n) = 0.5 ± 0.4 ± 0.1 (< 1.2), 
B{B" riTr") = 0.9±0.4±0.1 (< 1.5), B(B" r/Tr") = 0.9±0.4±0.1 (< 1.5), and S(B° ^ W7r°) = 
0.07 ± 0.26 ± 0.02 (< 0.5). The first error quoted is statistical and the second systematic. For the 
?7P+ mode, we measure the charge asymmetry Ach{B~^ —> rjp^) = 0.13 ± 0.11 ± 0.02. 
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PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



Measurements of charmless B decays are now rou- 
tinely used to test the accuracy of theoretical predic- 
tions based on, for example, QCD factorization [l|, 3|, 
flavor SU(3) symmetry [1 i, [| , perturbative QCD [6|, 
or soft coUinear effective theory [7| . We present measure- 
ments of the branching fraction and charge asymmetry 
for the decay rjp'^ (charge conjugate reactions 

are implied throughout this paper), superseding our pre- 
vious result that found a 4.7cr signal for this decay [§| 
with a luminosity of about one-half that used in this pa- 
per. In addition we search for the decays i?" rj'rj, 
B° r]Tr°, B° tj'tt^, and B° LUTr°. None of these 
decays have been observed previously though limits have 
been reported by BaBAR Q, Belle and CLEO [ll|. 

In the Standard Model (SM) the dominant processes 
that contribute to these decays are described by tree am- 
plitudes and to a lesser extent penguin (loop) amplitudes. 
For i?" rj'-K^ and B^ -qn^ the color-suppressed tree 
diagram is suppressed by approximate cancellation be- 
tween the amplitudes for the 7r° and for the isoscalar 
meson that contains the spectator quark. The approxi- 
mate ranges of expectations [H 0, y, 0, H, @, 0] for the 
branching fraction are ^10x10"^ for i?+ — > ?7/9+, 0.3- 
2x10-*^ for 5° ?/r/, 0.2-1.0x10-^ for t^('\° , 

and -0.1x10"^ for B° ^ ujn° . Direct CP violation 
could be detected as a charge asymmetry, defined as 
Ach = (r^^r^)/(r^+r^), where the superscript on the 
width r corresponds to the sign of the B^ meson; Ach 
for i?+ rjp^ is expected to be small since the decay is 
dominated by a single amplitude. 

These B'^ decays are also of interest in constraining 
the expected value of the time-dependent CP-violation 
asymmetryparameter 5/ in the B decay with final state 
/ = V'K^ BE El- The leading-order SM calculation 



S 



Jlii>K% 



where the latter has 



gives the equality S^ij^o^ 

been precisely measured [ij], and equals sin2/3 in the SM. 
The CP asymmetries in the charmless B decays are not 
only sensitive to contributions from new physics, but also 
to contamination from sub-leading SM amplitudes. Re- 
cent theoretical calculations of the size of the change in 
Srf'K° from these sub- leading amplitudes finds no more 
than 0.03 0, The most stringent constraint from 

data on such contamination uses SU(3) and the measured 



branching fractions of the decays B'^ 



B° 

gested [16|, ll7| that S° 



■q'-q, B° 



rjTT 



1% |13[ . Recently it has also been sug- 
j'tt'^ and B^ rjir'^ can be 
used to constrain the contribution from isospin-breaking 
effects on the value of sin2a in i? ^ tt+tt^ decays. 

The results presented here are based on data collected 
with the BaBAR detector 18| at the PEP-II asymmetric 



Center. We recorded a data sample at the T(45) res- 
onance (center-of-mass energy y/s — 10.58 GeV) with 



an integrated luminosity of 418 fb , corresponding to 
(459 ± 5) X 10^ BB pairs. 

Charged particles from the e+e" interactions are de- 
tected and their momenta measured by a combination 
of five layers of double-sided silicon microstrip detectors 
and a 40-layer drift chamber, both operating in the 1.5 T 
magnetic field of a superconducting solenoid. Photons 
and electrons are identified with a CsI(Tl) electromag- 
netic calorimeter (EMC). Further charged particle iden- 
tification (PID) is provided by the average energy loss 
{dE/dx) in the tracking devices and by an internally re- 
fiecting ring imaging Cherenkov detector (DIRC) cover- 
ing the central region. 

We establish the event selection criteria with the aid of 
a detailed Monte Carlo (MC) simulation of the B produc- 
tion and decay sequences, and of the detector response 
[l^. These criteria are designed to retain signal events 
with high efficiency. Applied to the data, they result 
in a sample much larger than the expected signal, but 
with well characterized backgrounds. We extract the sig- 
nal yields from this sample with a maximum likelihood 
(ML) fit. 

The i3-daughter candidates are reconstructed through 
their decays tt" 77, 77 — > 77 (7777), r/ Tr+Tr^Tr" (jysTr), 
■q' ri^^n+TT- (77^,^^), 77' p°7 (77' ), u tt+tt-tt^, 
— > TT+TT^ and /3+ 7r+7r°. Table [J lists the require- 
ments on the invariant masses of these particles' final 
states. All requirements are kept loose (> 3a) for later 
fitting except for the 7r° invariant mass which is not in- 
cluded in the fits. Secondary charged pions in 77, 77', and 
Lu candidates are rejected if classified as protons, kaons, 
or electrons by a combination of their DIRC, dE/dx, and 
EMC PID signatures. 



TABLE I: Selection requirements on the invariant masses of 
resonances and the laboratory energies of photons from their 
decay. 

State Invariant mass (MeV) E{'y) (MeV) 

Prompt 7r° 120 < m(77) < 150 > 50 

q-yj 490 < m(77) < 600 > 100 

qs^ 520 < m(7r+7r-7r'') < 570 > 30 

q'^^^ 910 < m{Tv+TT-q) < 1000 > 100 

q'p^ 910 < m(7r+7r-7) < 1000 > 200 

Lo 735 < 7n(7r+7r-7r'') < 825 > 30 

p° 510 < m(7r+7r~) < 1000 — 

p+ 470 < r7i(7r+7r°) < 1070 > 30 

We reconstruct the _B-meson candidate by combining 
the four-momenta of a pair of daughter mesons, with a 
vertex constraint if the ultimate final state includes at 
least two charged particles. Since the natural widths of 
the 77, 77', and tt*^ mesons are much smaller than the res- 
olution, we also constrain their masses to nominal val- 
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ues [20| in the fit of the B candidate. From the kine- 
matics of the T{AS) decay we determine the energy- 
substituted mass Toes — \J \^ ~ Pb ^'^'^ the energy dif- 
ference Ai? = Eb~\\/s, where {Eb,Pb) is the i3- meson 
4-momentum vector, and all values are expressed in the 
T{4:S) rest frame. The resolution in ttt-es is 3.0 MeV and 
in AE is 24-50 MeV, depending on the decay mode. We 
require 5.25 GeV < toes < 5.29 GeV and \AE\ < 0.25 
GeV (< 0.2 GeV for B° r]'r] and B+ r]p+). 

Backgrounds arise primarily from random combina- 
tions of particles in continuum e+e" qq events {q = 
u,d,s,c). We reduce these with requirements on the an- 
gle 6't between the thrust axis of the B candidate's decay 
products in the T(4S') rest frame and the thrust axis of 
the rest of the charged tracks and neutral calorimeter 
clusters in the event. The distribution is sharply peaked 
near | cos^tI = 1 for qq jet pairs and nearly uniform for 
B-meson decays. We require | cos^tI < 0.7-0.9 depend- 
ing on the decay mode. 

In the ML fit we discriminate against qq background 
with a Fisher discriminant !F that combines five variables 
[2H : the polar angles, with respect to the beam axis in 
the T{AS) rest frame, of the B candidate momentum and 
of the B thrust axis; the flavor tagging category [i^]; and 
the zeroth and second angular moments of the energy 
flow, excluding the B candidate, about the B thrust axis. 
It provides about one standard deviation of separation 
between B decay events and combinatorial background. 

We also impose restrictions on decay angles to exclude 
the most asymmetric decays where soft-particle back- 
grounds concentrate and the acceptance changes rapidly. 
We define the decay angle 9^^^ and its cosine Tik for a 
meson k as the angle between the momenta of a daughter 
particle and the meson's parent, measured in the meson's 
rest frame. We require for the 77^^ decays \HpO \ < 0.9 and 

for rj^'^TT^ iTi-n" I < 0.95. For i?° v'pjV-yj suppress the 
background B K*^ by requiring |7i^| < 0.86. These 
distributions are uniform for signal except for TipO which 
has a 1 — Hp distribution. 

For the 'qp'^ decay, we define 0^^^, as the angle 

between the 7r° and the negative of the B momentum in 
the p rest frame. We require —0.75 < Hp+ < 0.95. For 
the i?° wtt" decay, \'Huj\ is defined as the cosine of 
the angle between the normal to the lo decay plane (the 
plane of the three pions in the oj rest frame) and the flight 
direction of the B, measured in the lj rest frame. Both 
of these quantities have a Ti^ distribution for signal. 

The average number of candidates found per selected 
event is in the range 1.06 to 1.47, depending on the final 
state. We choose the candidate with the largest proba- 
bility for the fit to the B decay tree. 

We obtain yields for each channel from a ML fit with 
the input observables AE, toes, J^, ™fc, fc = 1,2 (the 
daughter invariant mass spectrum of the 77, 77', w, or p+ 
candidate), and Tik the helicity of the uo or p+ candidate. 



The selected sample sizes are given in the second column 
of Table |TT1 Besides any signal events, the samples con- 
tain combinatorial background from qq (dominant) and 
BB with b c, and a component from other charm- 
less BB modes that we estimate from the simulation to 
be no more than two percent of the sample. The latter 
events have ultimate final states different from the sig- 
nal, but with similar kinematics so that broad peaks near 
those of the signal appear in some observables, requiring 
a separate component in the probability density function 
(PDF). 

The likelihood function is 



N 



£ = exp(-5]y,)nE^^ 



V, imEs')V, {AE^)V, {r)V, (toI) V, {m^),V, (K,p+ ) 



(1) 



where N is the number of events in the sample, and, for 
each of the three components j, Yj is the yield of events 
and Vj{x'^) the PDF for observable x in event i. For the 
mode B'-^ rj'^^^ri^T^ we found no need for the charm- 
less BB background component. For the B^ Vy-iP^ 
and B^ — > wtt" decays we split the charmless BB PDF 
into components made from backgrounds with and with- 
out a p+. The factored form of the PDF indicated in 
Eq.[l]is appropriate since correlations among observables 
measured in the data are small. Distortions of the sig- 
nal yields caused by this approximation are measured in 
simulation and included in the bias corrections and sys- 
tematic errors discussed below. 

We determine the PDFs for the signal and charmless 
BB background components from fits to MC simulated 
events. Large control samples of B decays to charmed 
final states of similar topology [B^ £''^(i^"'"7r~7r°)7r+ 
and 5+ D'^{K'^Tr~Tr°)p'^] are used to verify the sim- 
ulated resolutions in AE and toes. Where the control 
data samples reveal small differences from MC, we shift 
or scale the resolution used in the ML fits. We develop 
PDFs for the combinatorial background with fits to the 
data from which the signal region (5.27 GeV < toes < 
5.29 GeV and |A£;| < 0.1 GeV) has been excluded. 

We use the following functional forms for the PDFs: 
sum of two Gaussians for 7^sig(™Es), '^sig bb(^-^)' ^^'^ 
the sharper structures in 7^bb(toes) and Vj{mk); lin- 
ear or quadratic dependences for combinatorial compo- 
nents of ^-^{mk) and for Vqq{AE); quadratic func- 
tions for VjilHuj]) and 'Pj{TLp+); and a Gaussian of dif- 
ferent widths below and above the peak, plus a broad 
Gaussian, for Vj{T). The qq background in toes is de- 
scribed by the function xVl — x'^ exp [—^(1 — x"^)], with 
X = 2mEs/-\/s and parameter ^. These are discussed in 
more detail in Ref. [23| and can be seen in Fig. [T] for the 
_B+ rjp'^ decay. 

We allow the parameters most important for the de- 
termination of the combinatorial background PDFs to 
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TABLE 11: Number of events A'' in the sample, fitted signal yield Ys in events (ev.) with statistical error, measured bias, 
detection efBciency e, daughter branching fraction product (H'^Oi ''■nd measured branching fraction B with statistical error 



for each decay chain, and the measured charge asymmetry Ach for the decay B 



VP 



For the combined measurements 



the significance S (with systematic uncertainties included), 
parentheses the 90% C.L. upper limits. 



branching fraction with statistical and systematic error, and in 



1\ /T ^ J ^ 

Mode 


AT { \ 

iV (ev.j 


Ys (ev.) 


Bias (ev.) 


e (%) (%) S (a) 


B (10-^) 




\ 


T7P+ 












y.u 


9.9 ± 1.2 ± 0.8 


0.13 ± 0.11 ± 0.02 




104609 


326tl* 


17 ±9 


ID. / 


oy .4 




1 n -1- 1 A 




0.07 ±0.12 




47918 


123t^^ 


13 ±7 


11.7 


22.6 




9.1 ± 2.2 




0.28 ±0.21 


■n'v 












1.4 


0.5 ± 0.4 ± 0.1 


«l-2) 






2191 


8.8t^;t 


0.9 ± 0.5 


25.6 


6.9 












896 


o 9+5.1 
o.^_4.1 


0.2 ± 0.2 


16.7 


4.0 




l.U_]^ 4 








39723 


n 7+12.2 
^■'-8.6 


0.0 ±0.5 


25.6 


11.6 












20672 




2.0 ± 1.0 


18.2 


6.7 










Tin" 












2.2 


0.9 ± 0.4 ± 0.1 


«1.5) 






9085 


18.61??:? 


4.4 ±2.2 


20.5 


39.4 




0.4 ±0.6 








4030 


23.3lli:i 


1.7 ±0.9 


17.3 


22.6 




1 3+°'^ 


















3.1 


0.9 ± 0.4 ± 0.1 


«1.5) 






3784 


20.6t«s:^ 


1.8 ±0.9 


22.5 


17.5 




1 1+0.5 








19789 


19 9+18.4 
1^.^^16.3 


2.7 ± 1.4 


18.9 


29.5 












39822 


9 ^+19.9 
^•^-16. 8 


0.5 ±0.5 


18.4 


89.1 


0.3 


0.07 ± 0.26 ± 0.02 


«0.5) 





vary in the fit, along with the yields for aU components. 
SpecificaUy, the free background parameters are most or 
all of the following, depending on the decay mode: ^ 
for TTiES) linear and quadratic coefficients for AE, area 
and slope of the combinatorial component for m^, and 
the mean, width, and width difference parameters for J-. 
Results for the signal yields are presented in the third 
column of Table [Til for each sample. 

We test and calibrate the fitting procedure by apply- 
ing it to ensembles of simulated experiments composed 
of qq events drawn from the PDF, into which we have 
embedded the expected number of signal and charmless 
BB background events randomly extracted from the fully 
simulated MC samples. We find biases of 0-17 events, 
somewhat dependent on the signal yield. The bias values 
obtained for simulations that reproduce the yields found 
in the data are given in the fourth column of Table [TTl 
Figure [T] shows PDFs and projections of subsamples of 
the data, enriched with a threshold requirement on the 
signal likelihood (computed without the variable plotted) 
for the 5+ rjp'^ decay. Figure [2] shows projections for 
the other four modes. 

We determine the reconstruction efhciencies, given in 
Table [III a-s the ratio of reconstructed and accepted 
events in simulation to the number generated. We com- 
pute the branching fraction for each channel by subtract- 
ing the fit bias from the measured yield, and dividing the 
result by the efficiency and the number of produced BB 
pairs [23] . We assume that the branching fractions of the 
T(4S') to B+B- and B°B° are each equal to 50%. Table 
im gives the numbers pertinent to these computations. 

We combine results where we have multiple de- 



cay channels by adding for each channel the function 
~2\ii{[C{B)/C{Bo)] (g) G{B;0,a')}, where Bq is the cen- 
tral value from the fit, a' is the part of the systematic 
uncertainty uncorrelated with other channels, and 
denotes convolution with a Gaussian function; the part 
of the systematic uncertainty common to all channels is 
then added in quadrature. We give the resulting final 
branching fractions for each mode in Table |TT] with the 
significance, taken as the square root of the difference 
between the value of —2\n£{B) (with only additive sys- 
tematic uncertainties included) for zero signal and the 
value at its minimum. The 90% confidence level (C.L.) 
upper limits are taken to be the branching fraction below 
which lies 90% of the total of the likelihood integral in 
the positive branching fraction region. 

The systematic uncertainties on the branching frac- 
tions arising from lack of knowledge of the PDFs have 
been included in part in the statistical error since most 
background parameters are free in the fit. For the signal, 
the uncertainties in PDF parameters are estimated from 
the consistency of fits to MC and data in control modes. 
Varying the signal-PDF parameters within these errors, 
we estimate yield uncertainties of 0.3-7 events, depending 
on the decay mode. The uncertainty from fit bias (Ta- 
ble |lT| includes the statistical uncertainty from the sim- 
ulated experiments added in quadrature with one-half of 
the fit-bias correction. We estimate the uncertainty from 
modeling the charmless BB backgrounds by accounting 
for the uncertainties in the knowledge of their branching 
fractions. These additive errors are the largest system- 
atic errors for the modes with small signal yield (but not 
B+ r]p+). 
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FIG. 1: Signal-enhanced projections for rj^jp'^ (left) and -qsnP^ (right) for (a,g) rriES, (b,h) AE, (c,i) (d,j) m^, (e,k) rrip, 
(f,l) Tip. The total {qq plus BB) background fit function is shown as blue dot-dashed, signal as magenta dashed, and the total 
as a solid blue line. These plots are made with a minimum requirement on the likelihood that has an efficiency for signal of 
15-35% while reducing the background by between two and three orders of magnitude. 



Uncertainties in our knowledge of the efficiency, found 
from auxiliary studies, include 0.4% x Nt and 1.5% x N^, 
where Nt and are the number of tracks and photons, 
respectively, in the B candidate. The uncertainty in the 
total number of BB pairs in the data sample is 1.1%. 
Published data provide the uncertainties in the B- 
daughter product branching fractions (0.7-3.9%). The 
uncertainties in the efficiency from the event selection 
are about 0.5%. 

We observe the decay i?+ -qp^ with a significance 
of 9 standard deviations. The branching fraction (9.9 ± 
1.2 ±0.8) X 10-6 and charge asymmetry 0.13 ±0.11 ±0.02 
are in good agreement with the theoretical expectations. 
We do not find evidence for the other four decays though 
the sensitivity of these measurements is now comparable 
to the range of the theoretical estimates. 
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